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Based on the analytic expression of free energy for infinitely long Ising strip with finite width joined 
antiperiodically on a variety of planar lattices, we show the existence of first-order phase transition 
at the critical point of Ising transition. The critical gaps of the transition are also calculated 
■ analytically by measuring the discontinuities in the internal energy and the specific heat. 
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I. INTRODUCTION 

In the study of critical phenomena, the investigations about finite-size effect have been very active0. These 
introduce more stringent test on the universality by including universal critical amplitudes and amplitude relations 0- 
The finite-size effect is also important in the comparison between theory and experiment. For example, it relies on the 
I . knowledge of finite-size effect to extract the bulk properties from the experimental results. In the study of finite-size 
effect, a semi-infinite system for which the effect can be singled out is always of fundamental importance. Thus, the 
-t— > system of an infinitely long Ising cylinder with finite circumference has been studied for some times To 
be more specific, we define the system as infinite extension along the x-axis and finite extension L along the j/-axis. 
To form an infinitely long Ising cylinder, we can join the finite width periodically or antiperiodically. For periodically 
joined circumference (PJC, hereafter), the condition imposed on the coupling of the boundary spin variable <J m ,£ is 
I ' cr»ri,L+i = Um,i with — oo < rn < oo. For antiperiodic joined circumference (AJC, hereafter), the condition becomes 
0m,i+i = —Cm,i- F° r the- system with PJC, the finite-size scaling behaviors of the shift of critical temperature S (L), 
defined as 5 (L) — |T max (L) — T c \ /T c with the specific-heat peak T max (L) and the bulk critical temperature T c , was 
known. The exact scaling form was shown to be 6 (L) = bhxL/ (L) 2 , where b is a function of coupling ratios and 
independent of the sizeQ- Recently, we studied the finite-size scaling of the distribution of partition function zeros 
of such a systemp|. The leading finite-size scaling behavior of the imaginary part of a zero labelled by j is given 
as Imzj (L) ~ L~ x l v with the correlation length exponent and the leading scaling behavior of the real part of 
the lowest zero (j = 1) can be written as |Rez! (L) — z c \ ~ 1"^"™. Here X zero is another critical exponent which is 
• closely related to the shift exponent A characterizing the scaling behavior of S (^1 flfl ITlL Il2l Il3| ] . Our results showed 
| that the system with PJC has v = 1 and X zero = 2. 

For the system with AJC, the influences of the finite-size effect on thermodynamic quantities are quite different 
from those for the system with PJC. It was shown that the critical point of two-dimensional Ising phase transition 
is one of the partition function zeros for the system with AJC, and the singularity associated with the zero at the 
critical point is a first-order phase transition for finite circumference 8] . As the length of circumference goes infinite, 
the effect of antiperiodic boundary condition becomes negligible, and the singularity changes to be a second-order 
phase transition. This is a very intriguing feature. It is known that antiperiodic boundary condition may enhance 
the appearence of the Block walls in spin configurations. In general, the existence of Block walls will increase the 
. free energy but suppress the energy fluctuations of the system. However, for an infinitely long Ising cylinder, based 
on the study of the cumulative distribution of partition function zeros we found that the existence of the zero mode 
of spin wave is solely responsible for the occurence of the non-analyticity at the critical point. We notice that the 
appearence of first-order phase transition does not contradict with the Mermin-Wigner theorem[lJ, 15]. Because 
the zero mode of spin wave always appear in the spin configurations, the system can not be viewed as an effective 
one-dimensional system with short-range Hamiltonian. First-order transitions are characterized by the first derivative 
of the free energy at the transition point. In this paper, we further analyze the nature of the first-order transition 
quantitatively by directly measuring the discontinuities in the internal energy and the specific heat for the systems 
with AJC on planar square, triangular, and hexagonal lattices with anisotropic couplings. 

This paper is organized as follows. In Sec. II, the analytic expressions of free energy density for an infinitely long 
Ising cylinder with PJC and AJC on a variety of planar lattices are given and discussed. In Sec. Ill, we perform explicit 
calculations for the internal energy and the specific heat. Our results show that there exists discontinuities locating 
at the critical point of two-dimensional Ising transition for AJC. Then, in Sec. IV we measure the discontinuities 
appearing in the internal energy and the specific heat explicitly. Finally, we summarize the results in Sec. V. 
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II. FREE ENERGY 



Consider an infinitely long strip with finite width L on plane triangular or hexagonal lattices. The two infinitely 
extended sides may be joined periodically or antiperiodically. Then, Ising model with ferromagnetic couplings between 
nearest neighbors is defined on such infinitely long cylinder. The corresponding form for the free energies per site per 
fcfiT was shown to be 
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with s = 1 for triangular lattice and s = 2 for hexagonal lattice. This result can be obtained either by direct 
derivation^ or by taking one side to infinite in the results of finite lattices with appropriate boundary conditions JtI 1 1 6j| . 

In the above expression, R is some lattice-dependent function of temperature and irrelevant to the physical singu- 
larity of the theory, and A^s (fi = 0, 1, 2, and 3) can be expressed in terms of particular functionals as 
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for the hexagonal lattice. Here the tj's are given in terms of the hyperbolic tangential functions as 
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with 77,- = Jj/fc^T, and Ji is the coupling constant given in the i'th direction as specified in Fig. 1. Note that the 
expression of triangular lattice with t 3 — gives the result of rectangular lattice. 

The A value of Eq. (1) is one-half and zero for PJC and AJC respectively. Thus, the exact forms of the free energies 
subject to the two types of boundary conditions differ only in the possible (p + A)-values. To simplify the notation, 
we define p = p + A and change the sum of Eq. (1) to be over the p- values. Then, p is integer ranging from to L — 1 
for AJC, and it is half integer from 1/2 to L — 1/2 for PJC. When the width of Ising strip L is extended to infinite, 
the sum of Eq. (1) becomes a simple integration and the two expressions of different boundary conditions reduce to 
the same result. In this limit, the critical temperature of the Ising transition is determined by the condition|l7|. 



A Q -A 1 -A 2 -A 3 = 0. 



The form of the free energy density of Eq. (1) can be rewritten as 
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where 



7 M = I 2^ In [l + G(p,fj)cos(^-V(p, »?))], (17) 
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nl , Al+ 2A 2 A 3 cos(2irp/L)} 1/2 

G{P ^ = Ao-Mcos^p/L) ' (18) 



and the angle ijj(p, rj) is defined by 



tan ip (p, rj) = — . 19 

A 2 + ^3 cos(2irp/L) 



We notice that the integration of Eq. (17) may be completed to simplify the expression. However, to calculate 
the derivatives of free energy we have to be cautious of the continuity of the integrand to avoid ill-defined results. 
Particularly, if any sigularity appears in the integrand, the integration becomes improper. Therefore, for concreteness, 
in calculating the derivatives of free energy we first proceed with the derivatives and then perform the integration in 
Eq. (17). 

III. DERIVATIVES OF FREE ENERGY 

The dimensionless internal energy density e and specific heat c are defined as follows, 

* = 7T-. ( 2 °) 

c = -"at = W (21) 

Note that the measurement of Eq. (20) depends on the reference scale which is one of the coupling constants 
indexed by i. However, the expression of the specific heat is invariant under the change of the reference scale since 
^(^) 2 = /3 2 (a//3) 2 foran yJ . 

The calculation of Eq. (20) contains the derivative of Eq. (17) which gives 
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where we neglect the label of reference scale, and the validity for interchanging the derivative with the integration 
is, however, only ensured at the points where the integrand is continuous for any Using the technique of contour 
integration, we obtain the result of Eq. (22) as 

As indicated, the analytic form of dl/drj given by Eq. (23) is valid only for positive function 1 — G 2 (p, rj). To 
inspect the domain of the functional values, we plot 1 — G 2 (p, rf) versus 77 for some p-modes in Fig. 2 which indicate 
that 1 — G 2 (p, rj) is a positive and concave function of 77 for a given p-mode and the absolute minimum locates at the 
point p = and 77 = r/ c with 1 — G 2 (0 7 r] c ) = 0. Here i] c denotes the 77 value at the critical point of ising transition. 
These features are summarized in the following lemma and corollary: 

Lemma 1 For the ferromagnetic couplings, 1 — G 2 {p 1 rf) as a function of r) is bounded from below for a given p-mode. 
Proof. The function of interest can be rewritten as 

N (p, 77) 
[A - Ai cos(27rp/L)] 2 

where the numerator N (p, rj) is given as 

N (p, 77) = A\ [cos(27rp/L) - B (t?)] 2 + C (r,) , (25) 
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with 

Bfo) . WIM) , (26) 

C(,) = MS-^-^-^i^W. (27) 

Since the denominator of Eq. (24) is strictly positive, the sign of the function is completely determined by the 
numerator TV (p, 77). For the ferromagnetic couplings, we have B (rj) > 1, which implies N (p,r]) > N (0,77) for any 
given temperature 77. Furthermore, according to the definitions of the A^'s given by Eqs. (2)-(5) we have 

N (0, 77) = [(ao - ai - a 2 - a 3 ) (ao - ai + a 2 + as)] 2 > 0. (28) 
Therefore, TV (p, 77) is bound from below so is the value of 1 — G 2 (p, 77). ■ 

Corollary 2 For p = i/ie minimum of 1 — G 2 (p, 77) as a function of 77 locates at the bulk critical point 77 = 77 c wit/i 
zero m/we, and for p ^ the minimum is larger than zero. 

Proof. The value of N (0, 77) given by Eq. (28) vanishes only at the critical point. This is due to the facts as follows: 
The critical condition of Eq. (15) is equivalent to 

ao — ct\ — ct2 — Qlz = 0. (29) 

Moreover, for the ferromagnetic couplings we have ao > at > for k = 1, 2, and 3. This leads the expression, 
ao — ai + 012 + as, in Eq. (28) to be strictly positive. Thence, among the minima of N (p, 77) for various p-modes 
there exist an absolute minimum TV (0, 77J whose value is zero. In other words for p ^ we have 

TV(p,77)>/V(0,77)>/V(0,77 c ) = 0. (30) 

Hence, the domain of the functional value of 1 — G 2 (p, 77) is specified by the conditions as 

1 -G 2 (p, 77) >0 forp^O (31) 

and 

1-G 2 (0,77) > 1-G 2 (0,77 C ) -0, (32) 

and the proof is completed. ■ 

According to the corollary, the case of 1 — G 2 (p, 77) < never happens for the case of ferromagnetic couplings. 
However, the situation, 1 — G 2 (p, 77) = 0, docs occur at one point for which the expression of Eq. (23) is invalid. Thus, 
the result of Eq. (23) fails only at the critical point 77 = ?7 C for AJC, and it is valid at any temperature for the system 
with PJC. With the exclusion of the exceptional point from the result of Eq. (23), we can obtain the internal energy 
e of Eq. (20) as 



e =e s + 2^X> £ M> (33) 
p 

with 
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D eM = - / =y m l^G M . (35) 

Likewise, using the definition given by Eq. (21) we can perform the derivative of e to obatin the specific heat c as 

c = c s + ^]T£ c (p,77), (36) 
p 
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with 




For the case of periodic boundary condition, our numerical results from the measurement of specific heat indicate 
that the shift of critical temperature, S(L), does obey the scaling form y|, 5 (L) = bhxL/ (L) 2 with b — 0.135 for 
the rectangular lattice, 0.225 for the triangular lattice, and 0.166 for the hexagonal lattice when the case of isotropic 
couplings is considered. For the case of AJC, the situation, 1 — G 2 (p,r]) = 0, as shown in the corollary, occurs at one 
point, p = and rj = r\ c . Then, the left and right derivatives of the free energy are discontinuous at this location, and 
this leads to the finite jumps at the critical point. Thus, the first-order phase transition occurs at r\ c for the system 
with AJC. 



IV. MEASUREMENT OF DISCONTINUITY 



To speculate how the finite jumps appear at the critical point, we may consider, for instance, Eq. (35) for the 
internal energy. By direct calculation, one can show that the quantity, d [lnG(0,?7)] /drj, vanishes at the critical point 
r\ c . In this manner, although the denominator, yl — G 2 (0, rj c ), also vanishes, the value of D e (0,i]) evaluated in the 
limit of approaching r\ c can be finite. The existence of the jumps is then related to the sign change of the functions, 
d[lnG(0,r))]/dr) and i/I - G 2 (0,7y), across the critical point. While the function d [In G(0, i])] jdr\ appears to be 
continuous function with nonzero slope intersecting the real line at the critical point, its value essentially reverses sign 
across the critical point. On the other hand, the function yl — G 2 (0, r\) must have its minimum at the critical point, 
and thence it does not change sign across the critical point. Thus, as approaching the critical point the left and right 
limits of Eq. (35) do not coincide for the zeroth mode and form the jump. Because the zero mode is absent from 
PJC, the discontinuity only appears in the case of AJC. 

To measure the discontinuity appearing in the case of AJC, we introduce the critical gaps, say A £ and A c for the 
internal energy and the specific heat respectively, defined as 



A £ = lim e(T) - lim e(T) , 

T->T+ T->T C ~ 

A c = lim c(T) - lim c(T). 

T-»T C + T^Tr 



(39) 
(40) 



Thus, A £ is the amount of latent heat per site released from the first-order phase transition. Since only the zero mode 
contributes to the discontinuity, according to Eqs. (34) and (36) we have 
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Then, we may extract the size-dependence from the critical gaps by introducing the gap amplitudes, et £ and a c , as 



A f = 



sL 



and A r = 



sL 



(43) 



Note that in the above equation we have isolated the the proportionality factor s between total number of lattice 
sites N and the length L along the finite edge. Owing to the inversional proportionality to L, the critical gaps in fact 
vanish in the thermodynamic limit for which the non-analiticity changes to be the second-order phase transition.. 

The explicit forms for the gap amplitudes can be expressed as functions of the coupling ratios = Ji/Jj, and the 
results read 
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for triangular lattice, and 
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for hexagonal lattice. Here the subscript c indicates that the functions of temperature contained in the paretheses 
are evaluated at the critical point of the Ising transition determined by Eq. (15). Note that the reference scales for 
Eqs. (44) and (46) are employed as J\ and Ji respectively, according to our convention shown in Fig. 1, and the 
permutation of the other coupling strengths leave the results invariant in each case. Again, the results for square 
lattice can be obtained from Eqs. (44) and (45) by putting t$ and r%i to zero. 

The numerical results of the gap amplitudes for various ratios of coupling constants are summarized in Tables. I 
and II. Since the expressions of Eqs. (44)-(47) all tend to be positive definite, by means of Eqs. (39) and (40), the left 
limit for the value of the internal energy or the specific heat at the critical point is always greater than the right limit. 
For both the triangular and hexagonal lattices, the gaps of the internal energy increase when the coupling constants 
are strengthened, however for the specific heat the behaviors appear to be rather complicated. Here, we record solely 
some of the typical values subject to the cases of isotropic couplings: For the triangular lattice a £ and a c are 6.0 and 
0.905212 respectively, while reducing to the rectangular lattices the values are 4.0 and 1.098589. On the other hand, 
the gap amplitudes for the hexagonal lattice are subsequently 3.464102 and 1.734378 for the internal energy density 
and specific heat. 



V. SUMMARY 



Based on the exact study of infinitely long Ising cylinders with the length of circumference L, we show analytically 
that there exists first-order phase transition for the case of antiperiodically joined circumference. The transition point 
locates at the critical point of the Ising transition. This peculiar phenomena arises because of the existence of the 
zero mode of spin wave. The soft mode of spin waves causes the free energy to behave less smoothly around the 
critical temperature, and this leads to the finite jumps in the derivative quantities, such as the internal energy and 
specific heat of the system. The finite jumps, referred as the critical gaps, are found to be inversely proportional to 
the width of the system. We then extract the size-dependent factors from the critical gaps by introducing the gap 
amplitudes. The magnitudes of the gap amplitudes for the internal energy and specific heat are also given explicitly 
in terms of the ratios of the coupling strengths. In general, the gaps of the internal energy increase when the coupling 
constants are strengthened, however for the specific heat the behaviors appear to be rather complicated. In the limit 
of L — > oo, the nature of the non-analiticity becomes the conventional two-dimensional Ising transition. 
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Figurel.The planar (a) triangular and (b) hexagonal lattices with directional coupling constants, J\, J2, and J3 
Figure2.The values of 1 — G 2 (p, rf) as a function of 77 for a given p-mode specified in the curve. Here we take L = 
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TABLE I: A list of some values of the gap amplitudes and the critical temperatures for the triangular lattice where conventionally 
Ji is used as the reference scale and the ratios of coupling constants are given as r 2 i = J2/J1 and r 31 — J3/J1. 



V2i r 31 



a t {i)(r 2 i,r zl ) a c (r 2 i,r 31 ) 



j_ 

16 
i_ 
16 
J_ 
16 
J_ 
16 
J_ 
16 
J_ 
16 
J_ 
16 
J_ 
16 
J_ 
16 
1 



16 



i .9749669067 

I 1.120609893 

i 1.359594319 

i 1.745136904 

1 2.368135472 

2 3.387294494 
4 5.080921897 
8 7.939611600 
16 12.83539772 
I 1.248671356 
i 1.471030695 
i 1.845058443 

1 2.464588209 

2 3.490273808 
4 5.202643928 
8 8.096230369 
16 13.05035679 
i 1.674907167 
i 2.034628498 

1 2.650888515 

2 3.690116886 
4 5.438377279 
8 8.398210101 
16 13.46287486 
i 2.383122015 

1 3.001777420 

2 4.069256996 
4 5.884122781 
8 8.964879141 
16 14.23051189 

1 3.640956906 

2 4.766244030 
4 6.699628670 
8 9.989370850 
16 15.59947050 
2 6.003554839 
4 8.138513992 
8 11.76824556 
16 17.92975829 
4 10.60355406 
8 14.76046755 
16 21.75350911 
8 19.71670567 
16 27.92219046 
16 37.89016754 



2.980303124 
3.127565685 
3.357757709 
3.697719840 
4.172052393 
4.794576091 
5.561554227 
6.453530816 
7.444318063 
3.265311339 
3.492162591 
3.839911133 
4.335206961 
4.990854252 
5.800546669 
6.742319380 
7.788226528 
3.724368061 
4.095394305 
4.637804574 
5.363639576 
6.261540830 
7.304442724 
8.461116081 
4.521379708 
5.165044415 
6.039616654 
7.122162939 
8.372913203 
9.753480244 
6.000000000 
7.172866022 
8.639380759 
10.32255458 
12.16074113 
8.866690995 
11.08035144 
13.65416899 
16.43931998 
14.54302699 
18.81937049 
23.54209421 
25.86585276 
34.25059835 
48.49746442 



.9945387423 

.9938210612 

1.001236145 

1.025135541 

1.070033044 

1.124851910 

1.158588497 

1.136952941 

1.046706223 

.9867447193 

.9857854706 

1.002261462 

1.044751239 

1.104142168 

1.147521663 

1.136541280 

1.054645265 

.970816527 

.9707838264 

1.003032167 

1.065812648 

1.124690886 

1.133261620 

1.067584184 

.9433175594 

.9477520819 

1.002111403 

1.079055707 

1.119838839 

1.084320289 

.9052117204 

.9199187375 

.9984555121 

1.079677481 

1.094445451 

.864293479 

.8938292729 

.9930189587 

1.071852827 

.8300549508 

.8743124974 

.9878571688 

.8066370472 

.8619605402 

.7926740283 
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TABLE II: A list of some values of the gap amplitudes and the critical temperatures for the hexagonal lattice where conven- 
tionally J 2 is used as the reference scale and the ratios of coupling constants are given as r 12 = J1/J2 and r 32 = J3/ 'J 2 . 



ri2 r 32 



V 2 ' c 



ffl E (2)(ri2,r 32 ) a c (r 12 ,r 32 ) 



j_ 

16 
J_ 
16 
J_ 
16 
J_ 
16 
J_ 
16 
J_ 
16 
_L 
l(i 
J_ 
16 
i_ 
16 
1 



16 



.1418239113 

I .2051068932 

i .3088086098 

i .4686005983 

1 .6558979475 

2 .7724015781 
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